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For Doppler-broadened media operating under double-double electromagnetically induced trans¬ 
parency (EIT) conditions, we devise a scheme to control and reduce the probe-held group velocity 
at the center of the second transparency window. We derive numerical and approximate analytical 
solutions for the width of EIT windows and for the group velocities of the probe held at the two dis¬ 
tinct transparency windows, and we show that the group velocities of the probe held can be lowered 
by judiciously choosing the physical parameters of the system. Our modeling enables us to identify 
three signal-held strength regimes (with a signal-held strength always higher than the probe-held 
strength), quantified by the Rabi frequency, for slowing the probe held. These three regimes cor¬ 
respond to a weak signal held, with the probe-held group velocity and transparency window width 
both smaller for the second window compared to the hrst window, a medium-strength signal held, 
with a probe-held group velocity smaller in the second window than in the hrst window but with 
larger transparency-window width for the second window, and the strong signal held, with both 
group velocity and transparency window width larger for the second window. Our scheme exploits 
the fact that the second transparency window is sensitive to a temperature-controlled signal-held 
nonlinearity, whereas the hrst transparency window is insensitive to this nonlinearity. 

PACS numbers: 42.50.Gy, 42.50.Ex 


I. INTRODUCTION 

Electromagnetically induced transparency (EIT) is a 
phenomenon whereby a medium can be switched between 
states of transparency and opacity through the control¬ 
lable application of a weak probe field |lj, and applica¬ 
tions of EIT include slow light Sfl , optical switching Q, 
and optical quantum memory [f|. Double EIT extends 
the notion of EIT from controllable creation of trans¬ 
parency at a given frequency to creation of two trans¬ 
parency windows at two different frequencies [7j. 

Double-double EIT (DDEIT) extends EIT even fur¬ 
ther to the case that each of a signal and a probe field 
experience two transparency windows with the trans¬ 
parency window for the signal controlled by a coupling 
field and a probe field and the transparency window 
for the probe held controlled by the coupling held and 
the signal held @], DDEIT introduces the possibility of 
controlling propagation of and interaction between two 
bichromatic helds. 

One advantage of controlling light is slowing it down. 
Slow light is especially important for optical communica¬ 
tion and for quantum information processing. For optical 
communication, slow light enhances light-matter interac¬ 
tion times and thereby leads to an increase in nonlinear 
interactions Ml, In quantum computing, slow light 
enables storage of the quantum state of light for a suffi¬ 
ciently long time to enable quantum memory [©]. 
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Our main objective is to propose a mechanism to 
slow the probe held in the second transparency win¬ 
dow of Doppler-broadened DDEIT. For this slowing to 
be achieved, we need to balance two competing require¬ 
ments. One requirement for slowing the probe pulse in 
Doppler-broadened EIT is to reduce the driving held in¬ 
tensity because the group velocity is proportional the 
driving-held intensity. On the other hand, the driving- 
held intensity must be sufficiently large to circumvent 
inhomogeneous broadening |10|-[l4j. 

We derive an analytical expression to enable us to hnd 
a parameter regime where these competing requirements 
can simultaneously be satished. Our analytical technique 
is based on approximating the Maxwell-Boltzmann ve¬ 
locity distribution for atoms by Lorentzian distributions 
over the narrow but relevant domain of small atomic ve¬ 
locities [l2j . We hnd that the nonlinear interaction be¬ 
tween the probe and signal held maintain the width of 
the second window constant for high Doppler width. This 
result permits us to lower the intensity of the signal held 
further without losing the EIT transparency window and 
get lower probe held group velocity at the second window 
compared to the hrst window. We apply our scheme to 
the case of 8 ' Rb under realistic experimental conditions 
and show that group-velocity reduction of the probe held 
is feasible. 

We present our work in the following order. In Sec. El 
we introduce the optical density matrix element describ¬ 
ing the optical properties of the transition |1) —► |4) and 
the resultant atomic susceptibility, and we use these re¬ 
sults to calculate the widths of the Lorentzian-shaped 
transparency windows and the corresponding group ve- 



2 



FIG. 1: Four-level tripod electronic structure with high- 
energy state 1 4) and lower-energy levels 1 3), | 2 ), and 11 ) in or¬ 
der of decreasing energy. Transitions are driven by probe (p), 
coupling (c) and signal (s) fields with frequencies w x and 
detunings 5 X with x£ {p,c,s}. Dephasing rates are 70 ; for 
i€{2,3}. 

locities for the probe field. 

Doppler broadening due to temperature is incorpo¬ 
rated into the expression for susceptibility in Sec. IIIII 
We solve this susceptibility numerically and analytically. 


Our analytical solution is based on ignoring quadratic 
dependence of the probe-field Rabi frequency and em¬ 
ploying a Lorentzian approximation for a narrow band 
around the Gaussian Maxwell-Boltzmann distribution. 
This approximate expression enables an intuition about 
how to control group velocities’ reduction at the second 
window. In Sec. IIIID1 we present the procedure to re¬ 
duce the group velocity in the second window. Finally, 
we summarize in Sec. El 


II. STATIONARY-ATOM OPTICAL 
SUSCEPTIBILITY 

Consider the closed fti atomic model scheme depicted 
in Fig.|T]with electronic level |4) coupled to the lower lev¬ 
els |1), 12), and |3) by three coherent fields, namely the 
probe, coupling, and signal Rabi frequencies, f2 p , O c , S2 s , 
respectively. The |1) -O 12), |1) -f4 |3), and |2) F4 |3) 
transitions are dipole forbidden. The three fields are de¬ 
tuned from the electronic transition frequency u t3 be¬ 
tween states I*} and |j) by 

<5 P :=w 41 -w P , (1) 

Sc :=w 42 - W c , (2) 

5 s •—^43 (8) 

respectively. 

The analytical steady-state density-matrix element 
( pi 0 ) solution for a stationary atom, to first order in the 
probe-field Rabi frequency, can be approximated from 
the exact expression Q as 
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( 4 ) 


where 

S xy '■= S x - 5 y (5) 

is the two-photon detuning. Details concerning the 

derivation of Eq. fl3J appear in Appendix [Al 

We verified this expression numerically for weak signal 
and weaker probe Rabi frequencies, i.e., for the condition 

|L! c | 2 »|fi s | 2 »|Dp| 2 . (6) 

The decay rates in Eq. 0 are 

: = ( 7 ) 

and the coherence decay rate is 

r k i = 7k + 7i- ( 8 ) 


The dephasing rate between the forbidden transitions is 
not zero; therefore, 72 = 7^2 and 73 = 7 ^ 3 . 

In our system, we impose the equal-population condi¬ 
tion 

Pu ~ P33 ~ 0.5. (9) 

Condition @ makes the equations approximately solv¬ 
able analytically as the equations of motion for popula¬ 
tion (IA 10 I) are effectively decoupled from the equations of 
motion for coherence (IA 11 I) . Condition © is achieved by 
incoherent excitation from ground state | 1 ) to the excited 
state 14) with constant pumping rate. Thus, the diago¬ 
nal matrix elements are held constant by conditions ©■ 
See Appendix IB II for details concerning the atomic pop¬ 
ulation in our scheme. The optical linear susceptibility 
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for an atomic gas in three dimensions with M the atomic 
density and di 4 the dipole moment is 


X 


(i) 

p 


P14 M\du \ 2 

V = - T. - 

toft 


( 10 ) 


We can substitute Eq. (J4J) into the numerator for Xp 1 ^ in 
Eq. (fTOD . which is complicated so we express Xp 1 ^ as 



ir/ 

2 (B 1 + 2i B 2 ) 


( Ci + 2i C 2 \ 
l A 1 -2iA 2 ) 


( 11 ) 


with the terms A 12 , B 12 , and C\ 2 explained in the fol¬ 
lowing. 

To simplify Eq. 0 , we fix the value ^44 = 0. This is 
always true because the atoms are trapped to the dark 
state leaving level |4) unpopulated. The population of 
the other three levels depends on the Rabi frequency of 
the driving fields. See Appendix IB 11 for more details of 
the dark-state analysis and state populations. 

The variables in Eq. m are 
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( 12 ) 


We now have expressions for the steady-state solu¬ 
tion 0 and the corresponding susceptibilities for the 
probe field (Ell- 

Expression m is used to calculate and plot the sus¬ 
ceptibility, whose imaginary part is shown in Fig. [21(a), 
and whose real part is shown in Fig. db). This absorp¬ 
tion plot clearly displays the first probe window cen¬ 
tered at <5 p = 5 C and the second EIT window centered 
at (ip — dg Sq. 


A. Linewidth and Group Velocity 


The linewidth of each transparency window 1 £ {1,2} 
is given by the half-width at half-maximum (HWHM) "l,. 
The half-maximum values >7 are determined first by find¬ 
ing the maximum ft, maXi and the minimum h m ; ni values of 
window 1 as shown in Appendix IC 11 and then calculating 


(13) 


By solving 


MXp^] = >7 (14) 





5 P (MHz) 


5 P (MHz) 


FIG. 2: (a) Absorption and (b) dispersion as a function of the 
probe detuning 5 P , with numerical (dotted line line), analyti¬ 
cal (solid line), and approximate linear equation (dashed line 
line) for 74 = 18 MHz, 73 = 10 kHz, 72 = 40 kHz, f! c = 74 , 
n a = O. 374 , Sip = O.O 574 , 8 S = 9MHz, and 8 C = 0. 


for S pc and <5 ps separately, "li and ~\ 2 are determined re¬ 
spectively, with (5 pc = 1i and (5 ps = ~\ 2 . In Appendix IC II 
we also discuss the requirements 

|1A| 2 > 7274 (15) 

and 

|H s | 2 > 7374 (16) 

to overcome homogeneous broadening [l5[ at the first and 
second windows respectively, and to observe the presence 
of the EIT windows. 

If conditions m and Cl are satisfied, the HWHMs 
of the first and second window are 



respectively. Probe dispersion is shown in Fig.[2}b). For 
detuning (5 p chosen at the center of each window, disper¬ 
sion is zero or close to zero. 

For each of windows 1 and 2, group velocity is [j] 


v 


g 


—, n g = (w 0 - A p ) 

n,g 


dRetxp^] 

c«p 


(18) 


for n g the group index, d cen the detuning at the center of 
each window (1 and 2), and wq the transition frequency 
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between levels |1) and |4). Detuning S cen equals 5 C at the 
first window and equals S s at the second window. The 
partial derivative of the dispersion in the denominator is 
determined by 

= nm Re[xp 1} (^ P )] - RcLxp 1 ^ = <5 cen )] 

d5 p S p ^<5 ce n <^cen 

<5cen 

(19) 

Therefore, the partial derivative of the dispersion 


9Re[xp 1} ] 


dS D 


v |fi c 


(7274 + l^cl 2 ) 


( 20 ) 


at the center of the first window and 


dRe[Xp 1] ] 

dS p 


s B 


V l^ s | 2 

(7374+ |«s| 2 ) 2 


( 21 ) 


at the center of the second window. Equations. (EDI) 
and m yield the slope of the tangent line to points 
S p = Seen as shown in Fig. (2}b). 

In Fig. [2Kb), the group velocity is shown to be approx¬ 
imately constant in each of the two EIT windows, which 
can be seen by the straight-line tangents. The group ve¬ 
locity scales inversely with slope so the ratio of group 
velocities for each EIT window is the inverse of the ratio 
of the slopes for each window. From Eqs. m and m 
and from Fig. [2jb) , the group velocity at the first win¬ 
dow evidently exceeds the group velocity at the second 
window for the given parameters. 

Under conditions (fl5l) and (fl6l) . the group velocity re¬ 
duces to 


2c |B C | 2 
77 w 14 


at the first window and to 


v 


82 


2c[ai 2 

7 w34 


( 22 ) 


(23) 


at the second window. Hence, for stationary atoms, the 
group velocities in both windows are linearly propor¬ 
tional to the intensities of the respective driving fields. 


III. DOPPLER-BROADENED OPTICAL 
SUSCEPTIBILITY 

At non-zero temperature atoms move randomly due 
to thermal energy. Thermal atomic motion leads to a 
spreading of the absorbed frequency due to the Doppler 
effect, which broadens the optical line profile and is 
known as Doppler broadening. 

In this section we solve susceptibility numerically and 
also derive approximate analytical expressions as a func¬ 
tion of temperature. These results are used to find the 


widths of transparency windows and also group veloci¬ 
ties of the probe field in each of the two DDEIT win¬ 
dows. Our approximate analytical technique is based on 
approximating the Maxwell-Boltzmann velocity distribu¬ 
tion for atoms by Lorentzian distributions over the nar¬ 
row but relevant domain of small atomic velocities [l2|. 
This approximation is valid as large velocities are suffi¬ 
ciently detuned so as not to affect the optics. 

In our scheme, the electromagnetic field passes through 
a gas of atoms at temperature T. Each atom of 
mass m has a velocity v obeying the Gaussian Maxwell- 
Boltzmann distribution 



(24) 


with v the component of velocity v in the direction of the 
three co-propagating signal, probe, and coupling fields. 

One effect of moving atoms is detuning of resonant 
frequencies due to the Doppler shift, which results in a 
velocity-dependent probe-field susceptibility x P (n). For 
our Doppler-broadened system, the susceptibility is thus 
av erag ed over the entire velocity distribution according 
to dh] 


X P := / X P (v)f{v)dv. (25) 

The velocity-dependent expression for susceptibility is 
obtained from Eq. (HD by the replacement 

S x S x + —-, x e {p,c,s}, (26) 

c 

for 

{ W 14 = Wo, x = P, 

W24, X = c, (27) 

w 34 , x = s, 

the atomic frequencies and c the speed of light in vacuo. 

Our scheme relies on neglecting Doppler effect on two- 
photon detuning <5 xy ©, which is achieved for the co¬ 
propagating fields driving approximately equal transition 
frequencies: 


wo = W14 ~ W24 ~ W34. (28) 

This choice is commensurate with our case of a 8 ' Rb gas. 
For this atom, we assign 11), |2), and |3) to the 55 , i/ 2 
level with F = 1, tof = 0 and F = 2, niF = {—2,0} 
respectively. Level |4) corresponds to the 5 Pi/ 2 level 
with F = 2 and tof = — 1. Therefore, the quanti¬ 
ties {<5 xy } in Eqs. (fill) do not change under Doppler 
broadening. 

Integration of Eq. (P25l) corresponds to a convolution of 
Lorentzian \p with the Gaussian profile, which is known 
as the Voigt profile E2- The Voigt profile can be solved 
numerically but is hard to solve analytically eh, mm 

0. 
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The lack of an exact analytical solution inhibits finding 
a simple expression relating the group velocity or width 
of each EIT window to Doppler broadening. Instead, we 
approximate the Maxwell-Boltzmann distribution bv a 
Lorentzian function over a narrow velocity domain [121 ] 
to obtain an approximate analytical expression for the 
optical susceptibility. This approximation is valid inso¬ 
far as we are interested in the optical response near the 
spectral center. 


A. Lorentzian line-shape approximation 


L{u)/f(v) (a. u.) 

,/"\ 

/ \ 

/ / * \ 

yj 

2 u) 0 v Vln2 

-1.4 -0.7 0.7 1.4 ^ 

FIG. 3: Plot of Lorentzian function (dashed line) and Gaus¬ 
sian function (dotted line) vs. normalized atomic velocity. 


In this subsection, we determine an analytical ap¬ 
proximation to the optical susceptibility for a Doppler- 
broadened system. Our approximation uses a Lorentzian 
fit to the Maxwell-Boltzmann velocity distribution over 
a narrow range of velocity. We use this approximation to 
show that the first probe transparency window is inde¬ 
pendent of the signal-field Rabi frequency and the second 
transparency window is nonlinear in the signal-field Rabi 
frequency. Furthermore, we derive the connection be¬ 
tween the transparency window and the Doppler broad¬ 
ening width, which is directly dependent on the temper¬ 
ature. 

The Lorentzian line-shape function [ill ] 

l ry = 4 =—( 29 ) 

V c / 1 / 7 F pj /2 _|_ 

is a function of the atomic velocity with Wl is the 
HWHM of the Lorentzian profile. To see that the 
Lorentzian (l29l) approximates the Gaussian (l24l) well over 
a narrow domain, we first write both functions as Maclau- 
rin series. The Gaussian (l24l) is approximated by 


The two expansions m and (1321) are approximately 
equal under the conditions that 

~ v^ W g ’ <34) 

for 

_l«^i v ^«l. (35) 

cW G 

Combining Eqs. m and (l34l) yields the connection be¬ 
tween the Lorentzian linewidth and the temperature. 
These conditions are satisfied near the center of both 
function profiles as shown in Fig. (3[ where the higher- 
order terms of Eqs. (1301) and (1321) have insignificant in¬ 
fluence. 

Integration of Eq. (l25l) using L(v) instead of /( v) has 
two terms evaluated with the contour integral using the 
residue theorem. The final optical susceptibility, includ¬ 
ing the Doppler broadening effect, is 

X P (S P ) = h(S p ) + I 2 (Sp) (36) 


f (— 2.^ _ a;g(\/In2) 3 2 

\ c / y/nWc (Py/wW^ 
^{y/hPlf 5 _ 
c 4 VFW 4 

with 


with <5p the detuning ©. The terms on the right-hand 
side of Eq. (l36l) are 


(30) 


IT) ypK 

~2 B 1 + 2LB 2 + W L 


(37) 


and 


W G := 

c V m 

the HWHM of the Gaussian profile and 

/VW0\ _ 1 _ __w l 2 
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(31) 
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ir] yfn Ci + i C 2 

~ J Bi + W L + 2i B 2 A 1 -W l - 2iA 2 

_ Ci + iC 2 

Ai + Bi + 2i(i?2 — A 2 ) 

irjW'Ly/n 

X Wl + 4 A\ - A\ + 4iAi A 2 ' 


(38) 


(32) The HWHM h of the first transparency window, and 
the group velocity for this window, depend on h(S p ) but 
not on hiSp) over the domain of <5 P pertaining to the first 
window. In the case of the second transparency window 
for the probe field, both Ji(<J p ) and J 2 (<5 P ) are non negli- 
gible for calculating the HWHM I 2 and group velocity. 
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FIG. 4: Plots of Imfxp 1 ^] and Refxp^] vs probe detuning 5 P 
at different temperatures for 74 = 18MHz, 73 = 10 kHz, 
72 = 40kHz, H c = 74 , fi s = O. 3574 , S s = 9MHz, and 5 C = 0. 
(a), (c) and (e) are Im[x P ] and (b), (d), and (f) are Re[x P ]. We 
set T= (1, 10, 100) K for (a),(b), (c),(d) and (e),(f) respec¬ 
tively, which is equivalent to Wl = (34.8, 110, 348)MHz, re¬ 
spectively. The dotted line-line corresponds to the analytical 
solution using the Lorentzian line-shape function, whereas the 
dashed line line is the numerical solution using the Maxwell- 
Boltzmann distribution function. 


In Fig. [H we plot the imaginary and real terms of the 
susceptibility xV as a function of the probe-field detun¬ 
ing (5 p at various temperature values based on the av¬ 
erage susceptibility (1551) for the Maxwell-Boltzmann dis¬ 
tribution function f(v) and for the approximation using 
the Lorentzian function L(v). At Low temperatures, for 
which the broadening is low, there is a discrepancy be¬ 
tween the two functions. 

At higher temperatures, for which 


tion m and condition 

|F1 C | 2 > 72^L (40) 

are required to observe the first transparency window. 
These conditions m and m eliminate the homoge¬ 
neous broadening and reduce the effect of inhomoge¬ 
neous broadening, respectively. At a temperature for 
which the Doppler broadening satisfies condition (1551) . 
and Eqs. m-m are a valid approximation, satisfying 
condition (1401) certainly implies satisfying condition (1151) . 

As shown in Fig. [4j the width “I 2 of the second trans¬ 
parency window is not noticeably affected by varying the 
Doppler width Wl- The reason for the robustness of “I 2 
is that the nonlinear interaction in I 2 , but not in / 1 , pro¬ 
tects the second window from deleterious temperature 
effects. Therefore, the strong-signal-field condition is not 
required to overcome Doppler broadening damaging the 
second transparency window. In other words, condition 

|a| 2 > 73 Wl (41) 

is no longer mandatory to observe the second window. 

Condition (fl 6 l) is still required to eliminate the homo¬ 
geneous broadening for significant transparency at the 
second window; see Appendix 1C II for detailed mathe¬ 
matical proofs of the conditions required to observe the 
transparency windows. Furthermore, the relaxation of 
condition (14T1) leads to further reduction of group veloc¬ 
ity in Doppler-broadened media, which was limited by 
the Doppler width appearing in the right-hand side of 
condition (14T1) . 

The two terms 72 Wl in Eq. (l40l) and 73 Wl in Eq. (14T1) 
quantify the inhomogeneous broadening of the two EIT 
windows. In other words, the Doppler broadening alone 
is not the whole story; rather, the products 72 , 3 Wl in¬ 
corporating the rates 72 and 73 are the key quantities. In 
Sec. IIII Bl we derive the linewidth and the group velocity 
for which the requisite conditions (flbl) and (1161) for elim¬ 
inating homogeneous broadening, are always satisfied for 
both windows. 


B. Width of the transparency windows 


W L 2 > 74, (39) 

the numerical data agree with the analytical data near 
the center as seen by comparing the two plots. The plots 
differ at the tail, which describes far-off-resonant atoms 
whose contribution is negligible. This numerical result 
validates the Lorentzian approximation for condition (1391) 
near the center, which leads to a rather simple form of 
the inhomogeneously broadened susceptibility. 

Analyzing the numerical result reveals that condi¬ 


The EIT width in a three-level Doppler-broadened A 
system can be maintained by keeping the temperature of 
the system constant while changing the driving field E 3 - 
fl4| . Here, we follow a different approach by studying the 
dependence of the linewidth on temperature while fixing 
the intensity of the driving fields. The intensities of the 
driving fields are chosen such to eliminate the homoge¬ 
neous broadening. 

The HWHM of the first window for the Doppler- 
broadened system is equal to 
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FIG. 5: Numerical (dashed line) and analytical (dotted line) 
solutions of the HWHM (“1) for the (a) first and (b) second 
EIT transparency windows vs Doppler width Wl for 74 = 18 
MHz, 73 = 10 kHz, 72 = 40 kHz, Q c = 74 , = 0 . 35 7 4 , 

<5 S = 9MHz, and <5 C = 0. Inset: Numerical (dashed line) and 
analytical (dotted line) HWHM of the second EIT window 
evaluated for the gain term eliminated. 


The width decreases nonlinearly as the Doppler 
width Wl increases as shown in Fig. [5j The condition 
|H C | 2 O 72 Wl is valid for all Wl values in the figure. 
For a high-intensity coupling field (l40l) . the width of the 
first window reduces to 


|^ c | 2 

2v / Wl(2 74 + W l )' 


(43) 


The formula for HWHM can be further simplified if 
Wl » 74 , thereby yielding 


_ l^cl 2 

* 1 "tWl" * 


(44) 


This result is consistent with the previous result for a 
three-level A atom, subject to a high-intensity driving 
field, for which the linewidth is proportional to the in¬ 
tensity of the driving field and inversely proportional to 
the Doppler width [121 ] . 

The HWHM of the second window of the Doppler- 
broadened system has a more complicated form than for 
the first window: 


.=1 _ |D S | 2 / (74 + Wl) + Wl [>c 2 (74 + Wi) — \ 

2 V 4 Wl 74 2 [3 - ^ 2(74 + Wl)] + |D S | 2 (74 + W L ) 


( 45 ) 


where 


2 73 W L + D 2 
(74 + Wl)( 7 3 Wl + H 2 ) 

2 73 Wl(74 — Wl) + fl 2 ( 2 74 — Wl) | 

W L [ 73 W£ + (Wl - 74)0 2 ] J 1 j 



is the half-maximum value of Imy p of the second window. 
The dependence of the HWHM of the second window 
on Doppler width is shown in Fig. [5] The width of the 
second window slightly decreases as the Doppler width 
increases. 

For large Doppler broadening, Wl O' 74 , m 2 de¬ 
pends on the population difference pn — p 33 . As we set 
P 11 ~ P 33 ~ 0.5, ir 2 is always located at Imx P « 0, i. e., 
where absorption vanishes. Consequently, the width of 
the second window remains approximately constant with 
respect to Doppler width 


% 


lai 2 

2'f2\J + 2 |D S | 2 


(47) 


This independence Doppler broadening width response 
of the second window is due to the gain described by 
Im / 2 of Eq. (l36l) . Expression (E| reveals that further 
reduction of the group velocity can be achieved by re¬ 
ducing the intensity of the signal field without losing the 
transparency window due to Doppler broadening. 

The two EIT windows have the same width at the in¬ 
tercept point between the two curves as shown in Fig. [5] 
For all values of Doppler width the signal field has lower 
intensity than the coupling field. The inset to Fig. [5] 
shows how the second window would behave as a function 
of Wl if the nonlinear contribution / 2 were suppressed. 
This inset makes clear how important the optical non¬ 
linearity is for achieving quite different temperature sen¬ 
sitivities of the two transparency windows for the probe 
field. Mathematically, an effect of forcing / 2 = 0 is that 
the HWHM of the second transparency window is given 
by a modification of the HWHM of the first window (l42l) 
with the proviso that fl c is replaced by D s and y 2 is re¬ 
placed by 73 . 

For atoms copropagating with the probe field, the gain 
term suppresses the narrowing of the width results from 
Doppler broadening. Generalizing the choice of atomic 
propagation direction relative to the direction of the three 
driving fields would of course lead to different results [HI . 

In summary, Eq. (1151) is the full expression of the 
HWHM of the second transparency window and accounts 
for the nonlinear interaction between the probe and the 
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FIG. 6 : Plots of the numerical (dashed line) and analytical 
(dotted) results for the partial derivative of dispersion with 
respect to Doppler width for (a) the first window and (b) the 
second window for 74 = 18 MHz, 73 = 10 kHz, 72 = 40 kHz, 
fl c = I. 574 , fl s = O. 574 , (5 S = 13.5 MHz, and 8 C = 0. Inset: 
numerical (dashed line) and analytical (dotted line) results 
for the partial derivative of dispersion vs Doppler width at 
the second window for 1 2 = 0 . 


signal field. Its behavior is depicted in Fig. [5] and shows 
the insensitivity of the second transparency window on 
temperature, which is represented by width Wl- Con¬ 
trariwise, the first window is sensitive to Wl- 


C. Group velocities at the transparency windows 


From Sec. IIIIA1 we have approximate analytical ex¬ 
pressions for susceptibilities at the two transparency win¬ 
dows. In this subsection we determine the derivative of 
the susceptibility with respect to the detuning <5 P and 
use these partial derivatives of dispersion (fl9l) to calcu¬ 
late the group velocities for the probe field in each of the 
two transparency windows. The response of the partial 
derivative of dispersion with respect to Doppler broad¬ 
ening system is shown in Fig. [ 6 j In this figure, numerical 
calculations show constant group velocity at the first win¬ 
dow and a sharply increased group velocity at the second 
window. 

The analytical expression for the group velocity of the 
Doppler-broadened system is evaluated using Eq. Cl 
but with the Doppler-broadened susceptibility (HUH) re¬ 
placed the free Doppler-broadened susceptibility Xp 1 -* [||. 
The partial derivative of Re[y p ] at the center of the first 
window is 


dRe[y p ] = 

96p (72WL + |ttc| 2 ) 2 


( 48 ) 


FIG. 7: Plot of the numerical (dotted) and analytical (dashed 
line) results for the partial derivative of dispersion vs coupling 
field at the first window for Doppler width Wl = 409 MHz, 
74 = 18 MHz, 73 = 10 kHz, 72 = 40 kHz. 


and at the center of the second window is 

12 


<9Re[x P 


d5 n 


277v6F|f2 s | 74 


*■ (74-WL)(73WL + |a| 2 )^ 

477A/7F | n s | 2 74 

W L |H S | 4 


(49) 


For the first transparency window and for a strong cou¬ 
pling field (l40l) . the group velocity of the probe field at 
the center of the first window has the same group veloc¬ 
ity as for the Doppler-free case (f22j) . The negligibility of 
the Doppler broadening effect is due to the intensity of 
the coupling field being large, as can be explained from 
the analytical expression (l48jl . 

Figure [ 6 ] shows agreement between the analytical ex¬ 
pression (1351) and the full numerical result applicable for 
small Wl- This agreement diminishes slightly as Wl in¬ 
creases. Therefore, the Lorentzian function can be used 
to study the Doppler-broadened dispersion response of 
the A configuration comprising the three states |1), 12), 
and 14) provided that condition (l40l) is satisfied. 

Our analytical expression is reliable in practical pa¬ 
rameter regimes. This agreement between the analytical 
Lorentzian approximation and the full numerical result 
under condition (1401) is presented in Fig. [7] for varying 
coupling-field Rabi frequency. 

We establish reliability of our approximation by com¬ 
paring to an approximate Lorentzian expression derived 
for a A EIT system Q. In our notation, their result for 
group index is 


n e oc 


_ 74 |»c | 2 

72 (74 + Wl) + |D C | 2 


(50) 


with the relation between group index and derivative of 
dispersion (1351) given by Eq. m- We can neglect 74 
from (1501) according to the approximation (1151) . Although 
result (15U1) is derived for a A system and our result (1351) 
for a ft) system, both results pertain to an EIT window in 
























9 


a strong-coupling regime, and the two Lorentzian-based 
approximations agree. 

At the second window, the analytical calculation fits 
the numerical solution for all chosen Doppler widths in 
the figure. Eliminating / 2 , (1551) leads to an equation for 
group velocity at the center of the second window being 
similar equation to Eq. (PP51) but with f l c replaced by D s 
and 72 replaced by 73. Similar dependence on Doppler 
width is shown in the inset of Fig. [G] 

To achieve matched group velocity for the probe pulse 
propagating through the first and through the second 
window, a non-zero nonlinearity is required. The nonlin¬ 
earity / 2 is zero only if the condition P 44 = P 33 = 0 is 
met. This case for nonlinearity is depicted in the inset 
of Fig. [G] By fixing P 44 = P 33 = 0 we have the unwanted 
additional effect of violating condition (TTB1) and thereby 
destroying the second window. 

The intercept point between the two curves shown 
in Fig. [6] reveals the operating temperature for group- 
velocity matching. At temperatures exceeding the 
matched group-velocity case, the group velocity in the 
first window is lower than the group velocity for the sec¬ 
ond window and vice versa for temperatures lower than 
the condition for matched group velocity. 

In summary, we demonstrate three important points 
in this subsection. First, the Lorentzian approximation 
is a useful and valid approximation for studying the dis¬ 
persion response of the probe field as long as the condi¬ 
tions (HU1) and (TUI) for / 2 = 0 [Eq. (1551) ] hold. Second, the 
second term of Eq. (1551) modifies the optical dispersion at 
the second window, which leads to a capacity for group 
velocity control through manipulating the temperature. 
Finally, due to nonlinearity, a signal-field intensity less 
than the coupling-field intensity does not necessarily im¬ 
ply that the probe field has lower group velocity at the 
second EIT window than at the first window. 


D. Group-velocity reduction 

In the previous Sec. IIII Bl and llH Cl we have studied the 
behavior of the width and the group velocity for both EIT 
windows of the probe field in Doppler-broadening me¬ 
dia. We have shown that a high-intensity coupling field is 
required to overcome inhomogeneous broadening, which 
represents an obstacle for group-velocity reduction. The 
width of the second EIT window is independent of tem¬ 
perature, which means that the enhanced group-velocity 
reduction is superior to the case that would hold if the 
width did depend on temperature as temperature depen¬ 
dence could only worsen this effect. 

In this section, we derive two expressions that relate 
the signal-field Rabi frequency f2 S to the coupling-field 
Rabi fl c and Doppler width W\ J . Satisfying the first ex¬ 
pression guarantees that the probe field has the same 
group velocity in each transparency window. Satisfying 
the second expression guarantees the same HWHM for 
the two EIT windows. 
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FIG. 8 : Plot of the partial derivative of dispersion (dotted- 
dashed line) and HWHM (dashed line) for the second EIT 
window and HWHM (upper horizontal dotted line line) and 
partial derivative of dispersion (lower horizontal dotted line 
line) for the first EIT window vs normalized signal-field Rabi 
frequency with fl c = 74 , Wl = 700 MHz, 74 = 18 MHz, 
73 = 10 kHz, and 72 = 40 kHz. 


The relation between S2 s and fl c can be satisfied for 
a wide range of temperatures bounded above and below 
by the requirements for the analytical approximations to 
be valid according to Eqs. in and dm. We then use 
these two expressions to divide the signal-field intensity 
to three regimes: a low-strength regime where the group 
velocity and EIT width are lower than the first window, 
a high-strength regime where both group velocity and 
width of EIT window are greater than for the first win¬ 
dow, and a middle regime where the group velocity is 
lower and the width is higher than for the first window. 

In Fig. [51 we plot the HWHM and partial derivative 
of dispersion for the second EIT window using Eqs. urn 
and m, respectively. We also plot the HWHM and 
partial derivative of dispersion for the first EIT window. 
Intercepts between lines show which signal-field Rabi fre¬ 
quencies yield matched HWHM or group-velocity condi¬ 
tions. Matched HWHM occurs at H Sl and matched group 
velocity occurs at f l Sh with f 2 s , lower than f 2 Sh , and 1 
and h refer to lower and higher values, respectively. We 
can choose values of fl s to control which of the two win¬ 
dows has higher HWHM and group velocity. 

We exploit our analytical expressions for the HWHM 
and the group velocity at the center of each window to 
find the lower and higher boundary values of the signal 
field. Equating Eqs. (PHI) and (H71) for real values of fl c 
and f l s gives us the lower boundary value of fl s : 



Similarly equating Eqs. (1481) and (PHI) gives us the higher 
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boundary value of Q s : 


^s h = o\ o73Wl + 


3W L 


19 


2t4^\ 
^L 2 73 ) ' 


(52) 


Equations (1511) and (l52l) reveal which signal-field strength 
should be selected to achieve either matched width or 
matched group velocity, respectively. 

At certain Doppler width, the boundary values of fl s in 
Eqs. (l5Tj) and (l52l) can be tuned by varying the coupling- 
field strength fl c . Both the matched group velocity and 
the matched HWHM have lower value as is reduced. 

In summary our four-level atom optical system can 
be operated at the second window in three different 
regimes depending on the signal-field strength. In the 
low-strength regime, the second window has very low 
group velocity compared to the first window but also 
has a lower EIT width. However, we can operate in 
this regime for lower group velocity as long as the width 
is resolvable experimentally. Alternatively, in the high- 
strength regime, the second window has a higher group 
velocity than for the first window, which makes this high- 
strength regime less desirable for low group-velocity ex¬ 
periments. 


IV. CONCLUSION 

We have achieved our objective of showing that the 
second DDEIT window has advantages over first window 
with respect to obtaining an enhanced reduction of group 
velocity. The presence of a nonlinear interaction between 
the probe and signal fields in optical susceptibility plays 
a crucial role in enabling temperature-controlled modi¬ 
fication of the optical response. At the second window, 
this term signifies the ability to reduce the narrowing of 
width and thereby yields increases of the group-velocity 
as the Doppler width increases. The modified optical re¬ 
sponse due to nonlinear interaction permits observing the 
second window for low intense signal field and promises 
for more group velocities’ reduction in the second EIT 
window. 

By identifying the signal-field boundary values H S1 and 
H Sh , we are able to identify the regime of the signal-field 
strength values that could result in slower group velocity 
than for the first window. The low-strength regime is the 
best for realizing low group velocity, but the EIT window 
could be difficult to resolve. The middle-strength regime 
is more robust in that the second EIT window is expected 
to be resolvable and the group velocity is expected to be 
low. The high-strength regime is less interesting as the 
group velocity is relatively high. 

Our approximate analytical calculation succeeds in de¬ 
scribing the optical response of the Doppler-broadened 
four-level optical system and helps in analyzing the sys¬ 
tem in the presence or absence of the nonlinear interac¬ 
tion. These analytical calculations also provide us with 
intuition of how the width or group velocities change in 


a Doppler-broadened system. Importantly, our analyti¬ 
cal expression helps us to study the relation between the 
coupling and signal fields and to achieve matching of ei¬ 
ther widths or the group velocities of the two windows. 
These conditions are not intuitively clear otherwise, and 
hence would be difficult to discern using only numerical 
calculations. 
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Appendix A: Atom-field Hamiltonian and equations 
of motion 

The semiclassical Hamiltonian for the atom-field sys¬ 
tem depicted in Fig. |T], within the dipole and rotating- 
wave approximations, is 

H(t) =H 0 + H dl {t) (Al) 

with the free-atom Hamiltonian being 

4 

H 0 = a tJ := |«) (j \, (A2) 

1 = 1 

and the driving interaction Hamiltonian is 

Hdi(t) =^(n p e lulpt ai4 

+ n c e iWct <724 + fl s e i “ sf (734 + H.c.^, (A3) 

with H.c. denoting the Hermitian conjugate. In the in¬ 
teraction picture, with respect to the free-atom Hamil¬ 
tonian (IA2I) . the atom-field system Hamiltonian has the 
form 

V(t) =^(n p e -1,5pt <Ji4 

+ H c e ^°*<T24 + f2 s e ’^ si <734 + H.c.^. (A4) 

The Hamiltonian involves oscillatory terms at different 
optical frequencies. Therefore, our next step is to find a 
Hermitian operator to transform the interaction Hamil¬ 
tonian to a rotating frame in order to eliminate the time 
dependence. Thus, we construct the rotating-frame op¬ 
erator 

A = 3<5 p (fii + (2<5p + <5 c )cT22 + (2(5p + 4)033 + 2<5 p <744 (A5) 

and eliminate time dependence by the following rotating- 
frame transformation 

H'(t) = e iAt/h V{t)e~ iAt / h - A. (A 6 ) 
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The resultant Hamiltonian is given by 


tion density matrix elements: 



+ H c (J24 + H s (J 34 + H.c. 


(A7) 


for 


Hq :— Sp C &22 T + Spd '44 (A8) 

and S xy given in Eq. (|5|). 

The resultant Lindblad master equation is 3 


P 


4 

~ ft [P, {a ZJ pa JZ - cr n p - pcrjj) 

i<3 

4 

+ J2~Y~^ a nP a JJ - a nP - P a jj) ( A9 ) 

J =2 


for p the state in the rotating frame. The Lindblad mas¬ 
ter equation includes both spontaneous emission and de¬ 
phasing, where 7 ^ is the decay rate of state |j) —> |z) 
and 7 is the dephasing of state |z). 

By substituting Eq. (IA7I) into (IA9I) . we obtain 10 opti¬ 
cal Bloch equations. Six more optical Bloch equations are 
obtained from complex conjugates of the six off-diagonal 
density matrix expressions presented in the following. 
The 10 optical Bloch equations are 


Putt) = ^-^72 + Wpcj Pi 2 (t) 

~ 2 i~^*Pl4(t) + ^pP24(t)] , 

Pl3(t) = (^-^73 + i^psj Pl3(t) 

~ 2 I - ^sPl4(t) + ^pP43(t)] , 

pu(t) = (^-^74 + Wpj pu(t) 

+ 2 [^cPl2 + ^sPl3 + ^p (pll(t) — P44(t))\ , 
P23(t) = -T 32 — il5 sc ^ P23(t) 

- ^ [ttcP43(t) - tt*P24(t)} , 

P24{t) = 2^ 42 P 2i ^ 

~ 2 {—QpP2l(t) + fl c (p44(t) — P22(t)) — H S P23] , 

P 43 W = ^-^r 43 - i<5 s ^ P43(0 

+ 2 [ — ^*P23(t) + (p44(t) — P33{t)) 

-n; Pl3 (t)]. (Aio) 

Now we present the equations of motion for the popula- 


Pllit) =72lP22(i) +73iP33(i) T 741 / 244(^0 


2 ^pP4i(i) QpPuif) 

P22{t) = — 721^22 W + l32P33(t) + 742^44^) 

- \ [~ fi cP24(0 + ^cP42it)} , 

P33(t) = - 731/033^) - l32P33{t) + 743/044(0 

- \ i~ n *sP34{t) + O s p 4 3(*)] , 

P44(0 = - 74/044(0 - ^[^c/O24(0 - ^cP 42(0 + ^sP34(0 


- H*p 43 (0 + H p pi 4 (0 - Hpp 4 i(0]. 


(All) 


In summary, we reprise the master equation for a single 
tripod (fti) atom driven by three detuned fields and obtain 
the requisite equations of motion for the density-matrix 
elements to solve the dynamics. 


Appendix B: Dressed-state analysis and atomic 
Population 

Here, we analyze the population of the fh atom for 
various driving-field strengths and different detunings of 
the probe field. In Sec. EU we investigate the atomic 
population behavior using the dressed-state analysis and 
numerical calculations. We obtain general expressions 
for atomic populations in steady state for three cases of 
probe-field detuning. In Sec. IB 21 we introduce an ana¬ 
lytical expression for the case we studied in this paper, 
corresponding to the signal-field strength being greater 
than the probe-field strength. In Sec. IB 31 we study the 
Doppler effect on the atomic population based on a nu¬ 
merical analysis. 


1. General discussion of atomic population 

Finding a general analytical expression for the atomic 
population using Eqs. (IA10I) and (IA11I) is not feasible due 
to the difficulties of decoupling the equations of motion 
of coherence from those of the population. Therefore, 
we analyze the dynamics of atomic population depend¬ 
ing on the interpretation from dressed-state analysis and 
numerical calculation of the atomic population. 

General expressions for the eigenstates of the Hamil¬ 
tonian (IatI) are complicated. For simplicity, we choose 
to find the eigenstates for the three following tractable 
cases. 


(i. Spc — 0 


The first case <5 pc = 0, or, equivalently, <5 p = S c . We 
allow S s to assume any different value. In this case one 
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of the eigenvalues Ai = 0 corresponds to eigenstate 

n* ^ 

' ' Vl^cl 2 + |fip | 2 Vl^el 2 + |Op | 2 

(Bl) 

This eigenstate is a dark state as it does not contain 
a contribution from state |4) and is not coupled to 
state 14). This is obvious from studying the total dipole 
moment d 40 for a transition from state |^d) to the bare 
state 1 4). 

If the magnitudes of the coupling field and probe field 
are appropriately balanced, the negative sign in the su¬ 
perposition of |1) and |2) [Ea. dBlI) ]. which forms the state 
|^d)j causes the transition moment (V , d|^ 4 d| 4 ) to vanish. 
Therefore, if the atoms are formed in this state there is 
no possibility of excitation to |4), hence no absorption. 

For the case that the coupling field is much stronger 
than the probe field (fi c f 2 p ), state | 1 ) is almost equiv¬ 
alent to |^d). Thus, atoms decaying to state |1) are 
trapped in this state and remain there throughout the 
interaction. The atomic probability of being in state |1) 


Pii(t), 5Im[pi 4 (t)] 
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|n c r + l^pl 2 


(B2) 


and being in state | 2 ) is 


P 2 = |<2|p D )| 2 = 


|fi c f + l^pl 2 


(B3) 


We numerically solve the master equation and plot 
atomic populations in Fig. [9] After a time of order of 
the radiative lifetime, the atoms should be trapped in 
the dark state |^d); which we verify by comparing the 
populations in Fig. [9] with the calculated dark-state pop¬ 
ulations. Vanishing of the probe-field absorption Impi 4 
supports the claim that the atom has decayed into a dark 
state. Furthermore, state |4) does not become populated. 

For Fig. [U)a), the dark state is equivalent to state |1) 
whereas, for Fig. [9|b), it is a superposition of states |1) 
and 12). The atoms are pumped into the state by com¬ 
bined action of coupling signal and probe field and spon¬ 
taneous emission similar to optical pumping mechanism. 
At steady state, the distribution of atoms depends mainly 
on the magnitude of the driving fields following the rule 
of Eqs. (IB2I) and (IB3I) . 

The other three eigenstates are 


(b) 

FIG. 9: Populations of state |1) and |2) represented by pn 
(dotted line) and p 2 i (dashed line) respectively, and absorp¬ 
tion of the probe field represented by 5 Im[pi 4 ] (solid line) as a 
function of time t evaluated numerically by solving the master 
equation, (a) Coupling field is stronger than the probe field 
with fl c = 74 and fl p = O. 374 . (b) Coupling and probe fields 
have the same strength with Sl c = fl p = 74 . The system is ini¬ 
tially prepared with p^ = 1 and p£$ = P 33 = P 44 = 0. The 
chosen parameters are 74 = 18 MHz, 73 = 10 kHz, 72 = 40 
kHz, fi a = O. 374 , S s = 0.5fl c , <5 P = S c = 0. 


In conclusion, detuning plays an important role in the 
distribution of atoms; thus when 5 pc = 0 and after a 
time of the same order of atom relaxation time atoms 
are trapped to the dark state and their distribution in 
the bare state | 1 ) and | 2 ) depends on the magnitude of 
f l c and flp, although the atoms are not prepared in the 
dark state. 


b. bps — 0 


flp |1) + fl c |2) + 13) + 2Aj|4) 

IlM = , , * ” , 2 (B4) 

y/|S 2 p| 2 + |n c | + lA^-Ipal 2 +4|A l | 2 

with eigenvalues A,, (1 £ 2, 3,4) where each {A^ is a root 
of the eigenvalue equation 

4A 3 —4A 2 (<5 ps + <5 P ) + A(45 ps 5p - |H C | 2 - |fi p | 2 - |H S | 2 ) 
+ i5 ps (| fl c | + | V2 p | 2 ) = 0. (B5) 


Now, we study the case that the probe and signal fields 
are at the two-photon resonance with a |1) ft- |3) tran¬ 
sition; i.e., <5 ps = 0. We allow the coupling detuning S c 
to assume any value. In this case, the Hamiltonian (1X71) 
has an eigenvalue A' = 0 with eigenstate 

K)= , ~ Qt 11)+ , ^ 13) (B 6 ) 

j | n 8 | 2 + | n P | 2 a /| h s | 2 + | h p | 2 
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and eignvalues A' (jg 1,2,3) with eigenstates 
|1) + A ^ c c |2) + fl s 13) + 2A'|4) 

i^> = / , n , av ; o 2 

\ (A'l SpijZ + l^pl 2 + l^ s l” + 4 ^ 2 


(B7) 


where each A' is a root of the eigenvalue equation 

4A' 3 -4A' 2 (6 pc + S p ) + A'(4 S pc S p - |O c | 2 - |ff p | 2 - |O s | 2 ) 
+ <5 pc (|fl s | 2 + |fl p | 2 ) = 0. (B 8 ) 


Pii(t), 5Im[pi 4 (t)] 



The eigenstate \’ip' T) ) is also a dark state as it does not 
contain a contribution from state |4) and is not coupled 
to state 1 4). 

Atomic populations for states | 1 ) and |3) are calculated 
numerically and shown in Fig. [TIB At steady state the 
atoms are trapped in the dark state | ip' D ) as long as the 
coupling field is greater than or equal to the probe and 
to signal field. We claim that the atom is trapped in the 
dark state because, if instead the atom were in any one 
of the bright states of Eq. (lB7l) . the following phenomena 
would arise: 


Pii(t), 5Im[p 14 (t)] 



(i) We would expect to see some population in 
states |2) and |4) whereas, in Figs. IToT al and (b), 
the populations of states |1) and |3) add almost 
to one, hence making the combined population of 
states |2) and |4) nearly zero. 

(ii) For the case that Sl c fl s > fl p as shown in 
Fig. [He), if the system is in a bright state, then 
the population in state |3) will exceed the popula¬ 
tion in state 11 ), i.e., p 33 > pn, but the opposite is 
true: most of the population has been transferred 
to | 1 ). 

(iii) The absorption would not vanish for a bright state, 
but, in Figs. unw-nnrc), absorption vanishes; 
hence, the atoms are trapped in the dark state \if>' D ). 

At steady state the populations in | 1 ) and |3) are gov¬ 
erned by the signal- and probe-field Rabi frequencies ac¬ 
cording to 


p i = I( 1 IV’d )| 2 = 


l a 


|A2 S | + |fi F 


(B9) 


Pii(t), 5Im[pi 4 (t)] 
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(c) 

FIG. 10: Populations of level |1) and |3) represented by pn 
(dotted line) and P 33 (dashed line), respectively, as a function 
of time t evaluated numerically by solving the master equa¬ 
tion. (a) Signal- and probe-field strengths are of equal mag¬ 
nitude less than coupling field, with Q. s = Q. p = O .574 while 
fl c = 74. (b) Coupling-, signal- and probe-field strengths are 
of equal magnitude fi c = fi s = fi p = O. 3574 . (c) Signal- 

field strength is stronger than the probe field, with fl c =74, 
fi s = O.574, and ftp = O.I574. Other parameters are 74 = 18 
MHz, 73 = 10 kHz, 72 = 40 kHz, 5 a = <5 P = 0.5fl c , and <5 C = 0. 
Initial populations are p) ( f = 1 and = p^ = p^ = 0. 
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and 


^3 = l(3|^)| 2 = 


| Op 


l a 


I^pI 2 


(B10) 


respectively. 


C. dps — dpc — 0 

The last case pertains to the three detunings are 
equal Q, (d p = S c = S B ), which results in zero two-photon 


resonance. Two of the eigenstates are degenerate eigen¬ 
states with eigenvalues Ai = A 2 = 0 : 

I^di) - ~° s * |1> + ^ |3) 

y'lai +ia>l 2 viai 2 + ia>l 2 

- Tl c fl p |1) — (Op + fig) |2) + fl c fl s 13) 

llW = , ^=- 

J(iai 2 + iai 2 + iai 2 ) (ia ,| 2 + iai 2 ) 

(Bll) 
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These two states are dark as neither contains contribu¬ 
tion from state |4) nor involve transitions to state |4). 
However, the rest of the eigenstates retain a component 
of all the bare atomic states: 

= H p |l)+H c |2)+a|3)±2A ± |4) 

^lap + i^pP + iap + ^Ai) 2 

with 

A ± = \ (<5p ± ^ + o 2 + a + o 2 ) 

the corresponding eigenvalues. 

The steady-state atomic populations for the case 


(B12) 


(B13) 


Ai(t) 
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0.2 



t(s) 


0.0005 0.001 0.0015 

(a) 


t(s) 


<5 P = S a = S c = 0 (B14) 

are shown Fig. [TlJ In all cases the atomic population 
is distributed between states 11), |2), and |3) and ex¬ 
cludes |4). This exclusion suggests that, for all cases (a), 
(b), and (c), atoms are trapped in the dark state \ 1 pn 2 }, 
but we see now that this guess could be true for case (a) 
but not for cases (b) and (c). 

In Fig. HTT bl. we have O p , f l c ^ O s which means that, 
if the system is in dark state 12 / 702 ), then the population 
in |1) must be higher. However, we see that the popula¬ 
tion of state 13) is in fact higher and exhibits the opposite 
behavior to that shown in Fig. Etc). Thus, the system 
corresponding to Figs. fllTbl andfTITc) could be trapped 
in |'0 di)- However, the few populations in state |2) pre¬ 
vents us from making this conclusion as well. 

From this argument, we conclude that the system is 
not in a pure dark state, but relaxes into a mixture of 
the two dark states (IB11I) which is also a dark state given 
by 


PD = PDl \lpUl) (V’Dll + PD2 |V>D2) ('0D21 , (B15) 


where pui is the probability of being in state IV’di) and 
Py )2 is the probability of being in state \ 1 p-D 2 ) such that 

Pdi+PD2 = 1- (B16) 

The probability for state |1) being populated is 


A =(1 Ipd| 1) 

_ 1 
“|O p | 2 + |O s | 2 



p D2 |»p1 2 |»c| 2 A 

|Oc| + |o p | 2 + |o s | 2 y 

(B17) 


for populating state |2) is 


A = (2 |p D I 2) 


Pd 2 

Op + Hg 

,2 

^|O c | + |O p 2 + |O s | 2 j 

'( 

|Op| 2 + |O s | 2 ) 


(B18) 
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FIG. 11: Populations of levels 11), |2) and |3) represented 
by pn (dotted line), P 22 (dotted-dashed line), and P 33 (dashed 
line), respectively, as a function of time t evaluated numer¬ 
ically by solving the master equation. The conditions are 
(a) fi s = fl p = O .374 and fl c = 74 , (b) fl p ,fl c with 

fi c = fi p = 74 and fl s = O. 374 , and (c) tip with 

Q c — Q s = I 74 and fl p = O. 374 . The initial population is 
p® = 1 and pj °2 = P 33 = P 44 = 0. Other parameters are 
74 = 18 MHz, 73 = 10 kHz, 72 = 40 kHz, 5 P = S B = 5 C = 0. 
Insets (a), (b), and (c) are the absorptions represented by 
Impi 4 . 
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and for populating state |3) is 


P 3 =(3|pd|3) 

1 

- in p l 2 + iai 2 


^PDl |^p| 2 


+ 


PD2|a| 2 |»c| 2 \ 

iai + |iipi 2 + lap / 

(B19) 


The relation between P \, P 2 , and P 3 is 


Pi + P 2 + P 3 = 1 (B20) 


from TrpD = 1 . 

We use numerical calculation of the population in 
state 12) and Eq. (IB18I) to determine the value of p^ 2 - 
Once pd 2 is known, ppi is calculated from Eq. (IB 161) . 
The agreement between the numerical values of pn 
and P 33 and the calculated values of Pl and P 3 using 
Eqs. (IB 181) and (IB19I) . respectively, verify that the sys¬ 
tem is in a mixture of the two dark states (IB111) . 

We see that, for certain two-photon detunings, the sys¬ 
tem is eventually trapped in a dark state even if the atom 
has not been prepared initially (at t = 0 ) in a dark state. 
The atom is driven into the dark state by a combined 
action of coupling, signal, and probe fields and by spon¬ 
taneous emission similar to that occurring through an 
optical pumping mechanism. 

For stationary atoms (at low temperature), the steady- 
state atomic population depends on probe-field detun¬ 
ing due to dark-state dependence on probe-field detun¬ 
ing. Thus, changing the probe-field detuning modifies 
the steady-state population in each energy state, as long 
as the probe field has comparable strength to the sig¬ 
nal field even if both are quite weak compared to the 
coupling field strength as shown in Fig. fl^T a'). However, 
the dependence of the atomic population on probe-field 
detuning decreases as the probe-field strength become 
weaker than the signal-field strength; this feature is ev¬ 
ident by comparing Fig. UKa) with Fig. 112T b'). Almost 
all of the population is evidently trapped in the dark 
state |^d) when <5 pc = 0 and in dark state | V^d) when 
<5 ps = 0 . which corresponds to state | 1 ) when fl s f 2 p . 


2. Atomic population for probe-field strength 
weaker than signal-field strength 


In this section, we derive an analytical expression for 
atomic populations for the case studied in our paper cor¬ 
responding to fl c H s ^$> H p . The analytical expression 
can be found by solving Eqs. (IA10I) and (IA11I) restricted 
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FIG. 12: Numerically evaluated steady-state populations pn 
(dotted line), p 2 2 (dotted-dashed line), P 33 (dashed line), and 
P 44 (solid line) vs probe-field detuning <5 P evaluated (a),(b) at 
zero temperature and (c),(d) at 100 K. Parameter choices are 
(a) fl s = flp = O. 374 . (b) » 2 P , = O. 374 , 2 P = O.OI 74 , 
(c) n a = 2 P = O. 374 , and (d) 2 S 2 P , 2 S = O. 374 , 2 P = 

O.OI 74 . Other parameters are 74 = 18 MHz, 73 = 10 kHz, 
72 = 40 kHz, and 2 C = 74 , <5 S = 0.52 c , <5 C = 0. 


to the case that H p = 0 : 


P23(t) 

P 2 i(t) 


P43 {t) 


— ^32 — i^sc ) p23(t) — — f2 c P43(^)j 


--r 42 + Wc ) P2i{t) 


- 2 [^c (P44:(t) — P22(t)) - H s p 23 ] , 


-^r 43 - i<5 s ) p 4 3 (t) 


+ 2 [~^*P23(t) + ^s(P44:(t) — P3 3 (0)] ) 
Pllit) =l4lP4i(t) 

P22{t) =742P44 (t) - ^ [-0*p 24 (f) + H c p 42 (t)] , 

P33{t) =l43P44(t) - - [-fl*p 34 (t) + O s p 43 (f)] , 

1 =Pll(t) + P22(t) + P33(t) + p 44 (t). (B21) 


As we mentioned earlier the |1) 77 |2), |1) 77 |3), and 
|2) 77 |3) transitions are dipole-forbidden. Therefore, we 
restrict y 2 i = 731 = 7 32 = 0 in the above equation. 

With initial population described by 
Pii(0),P22(0),P3 3 (0), and p 44 (0), the atomic popu¬ 
lations for the four atomic bare states in the steady state 
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are 


P 11 

P 22 

P33 

P44 


la ZY 

741 YY + 743 (X-Y)’ 

741 ^^ 22 ( 0 ) + [Z(Y — 742 ) — X 7 43 )(l — pn(0)] 
74i2F + 74 3 (X + y) 

_ X 741 ^ 33 ( 0 ) 

741 ZY + 743 (X - Y) ’ 

Y \{Z + 743 )(1 ~ Pn(0)) + 74 iP 33 ( 0 )] 

74izy + 7 4 3 (x-y) 

(B2 

JaZY + ^iX-Y)' ' 


with 


X = 


\n s \ 2 \n c 


(rl 2 + 452 ) (r| 2 + 4<5 S 2 C ) 

r 43 + (4<5 S <5 SC + r 32 r 4 2 ) 


0 


.0 


43 


|o c | 2 r 3 2 

r++4+. 


■)+<( 


|n c | 2 ^c 

r++4<5 2 


- 6 , 


+- 


(i l sSi- i -) (4r 3 2«=-r 12 i c ) 


( v , |o c | 2 r 3 2 y , , ( |nd 2 5ac r y 

l i43+ + H4+^ _ds J J 


Y = 


|O c 


- 42 


r 2 

1 42 


- 4(5 2: 


(B23) 


and 


Z = 




l+l 2 r 32 A 

r| 2 +4 si) 


|£+| 2 r 32 ^\ 2 
rf 7+4Sf c ) 


+ 4 


( |n c | 2 f ac 

V r 32+4.5 2 c 



(B24) 


Equations (IB22I) tell us that, for all probe-field detun- 
ings, almost all the atomic population is in state | 1 ) with 
almost no population in state |3). This lack of population 
in 1 3) eliminates the effect of the nonlinear signal-probe 
interaction described by the second term of Eq. dU). As 
we require population in 1 3), we introduce an always-on 
incoherent pump at rate r to maintain population in |3). 

The equations of motion of the density matrix elements 
with the incoherent pumping are similar to those without 
incoherent pumping Eqs. (IA10I) and (IA11I) . differing only 
in the replacement 


74 -+74 + 2r, 73 -+ 73 + r, 

72 -+72 + r, r 34 -+ r 34 + r. (B25) 


The atomic population equations in the presence of the 
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FIG. 13: Steady-state populations pu (dotted line), P 22 (dot- 
dash), P 33 (dashed line), and P 44 (solid line) vs probe-field 
detuning 5 P for fi s fl p in the presence of incoherent pump¬ 
ing with constant rate r evaluated numerically by solving the 
master equation. The solutions correspond to (a) low tem¬ 
perature for r = 1MHz and (b) for two temperatures 100 and 
400 K with r = 0.01MHz. Other parameters are fl s = O. 374 , 
Q p = O.OI 74 , 741 = 742 = 6MHz, 743 = 12MHz 73 = 10 kHz, 
72 = 40 kHz, and fl c = 74 , S a = 13.5 MHz, 5 C = 0. 


incoherent pumping are modified as 


P 11 


p 22 

P33 

P44 


_ ZY (741 + r) _ 

Z (-r 742 + 4 rY + 741Y) + r7 43 (Y - X) ’ 
r(Z(Y- 142 )- 14 3 X) 

Z (—r7 42 + 4 rY + 741Y) + r7 43 (Y - X) ’ 
_ rY (Z + 743) _ 

Z (-T 742 + 4 rY + 741Y) + 4-743 (Y - X) ’ 

—- rZY - (B26) 

Z (-r 742 + 4 rY + 741Y) + r74 3 (Y - X) 


with X , Y, and Z defined in Eq. (IB23I) with replace¬ 
ment (IB25I) . 

The modified atomic population in the presence of in¬ 
coherent pumping for the case Q s O p is shown in 
Fig. 1131 a). The existence of incoherent pumping makes 
the population constant for all probe detunings. The 
value of the pumping rate r controls the population in 
each state. We use r = 1 MHz to populate states |1) 
and |3) with pu = p 33 = 0.44. At high temperature, 
when the Doppler effect plays a critical role in repopu¬ 
lating the states, |3) can be repopulated to a value of half 
using pump rate r as low as 1kHz as shown in Fig. 1131 b). 

In this section, we have presented an incoherent 
pumping procedure that maintains equal population be¬ 
tween |1) and 13) for a weak probe held. In the main 
body of the paper, we did not treat pumping; instead, 
we assumed equal population between | 1 ) and |3). As we 
use an incoherent pump, we should be concerned that 
coherence is affected, but we see here that dephasing 
due to incoherent pumping is negligible for reasonable 
parameters. Specifically, the extra dephasing of system 
from incoherent pumping is of the same order as 7^3 for 
Doppler-broadened media [20j |. 
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3. Temperature dependence of atomic population 

Increasing the temperature adds two more phenomena 
to the atom-field system which we incorporated into an 
extended quantum master equation. These two phenom¬ 
ena are thermal dissipation and Doppler broadening. Our 
examination of thermal dissipation shows that its effect is 
too weak to influence substantially either the coherence 
or the population. However, the second phenomenon of 
Doppler broadening modifies the coherence, as presented 
from Sec. EH onward for equal populations in levels | 1 ) 
and 1 3), is discussed in this subsection in the absence of 
this equal-population restriction. 

Finding an analytical expression for population in 
Doppler-broadening medium is difficult. Therefore, we 
perform numerical studies of atomic populations for var¬ 
ious temperatures. Now, we proceed to analyze the con¬ 
nection between atomic population and coherence. 

At high temperature, Doppler broadening reduces co¬ 
herence 0 , 0 ] and specifically directly reduces the co¬ 
herences pi 4 and P 34 that are established by the weak 
fields. Consequently the populations of the states |1) 
and 13) change according to the solution of Eq. (IA11I) . We 
see from the approximate analytical expression for opti¬ 
cal susceptibility (l36l) that increasing Doppler width Wl 
in Eq. (I5S1) is responsible for reducing coherence. 

Reduction of coherence and its consequent effects due 
to Doppler broadening and to increasing Doppler width 
are similar to the effects due to adding extra dephas¬ 
ing 7^41 between |1) and |4) and 7^43 between |3) and |4) 
plus increasing the dephasing 702 between |1) and |2). In 
Fig- El we present numerically evaluated atomic popu¬ 
lations for different probe field detunings d p at zero tem¬ 
perature accompanied by additional dephasing quantified 
by 7041 and by 7043 . We choose the dephasing 7041 and 
7043 in Fig. Q31 to be of the same order of magnitude of 
Wl- We observe that Figs. me), and I12f dl are similar 
to Figs. HiTa). andflflb). 


a. Us — Dp 

At <5 p = S c = 0, the population of atoms in state |3) 
increases as temperature increases, as shown Fig. me). 
Increasing the temperature from 0 to 100 K decreases p\\ 
from 1.0 to 0.7 while increasing p 33 from 0.0 to 0.2. This 
population changes is due to reduction of the coherence 
P14. The coherence P34 does not have an influence at 
<5 P c = 0 . 

The atom-field system is no longer in the pure dark 
state IV'd)- Increasing the temperature has the effect of 
displacing the system from the dark state to different 
state, where the absorption of the probe field increases. 
Thus, as the temperature of the system increases, probe- 
field absorption becomes very high, thereby potentially 
preventing the first EIT window from being observed. 

For <5 p = S s , the atomic population remains the same 
at different temperatures; i.e., pn = P 33 = 0.5. There¬ 
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FIG. 14 : Numerically evaluated steady-state populations pn 
(dotted line), P22 (dot-dash), P33 (dashed line), and P44 (solid 
line) vs probe-field detuning < 5 P at zero temperature and incor¬ 
porating dephasing 7041 = 7043 = 150 MHz and 702 = 0 . 8 MHz. 
The plots show (a) D s = D p = O.374 and (b) D s D p , 
D s = O.374, D p = O.OI74. Other parameters are 74 = 18 
MHz, 73 = 10 kHz, and D c = 74, 8 S = 0 . 5 D C , < 5 C = 0 . 


fore, the system remains trapped in the dark state |'i/'p) 
[Eq. (IB 6 I) ]. The population at that detuning is less 
sensitive to Doppler broadening. This insensitivity can 
be explained as resulting from higher-order nonlinear in¬ 
teractions between the signal and probe fields resulting 
from the coupling of P 14 to P 34 through the presence of 
the term pi 3 @. This coupling eliminates the effect of 
Doppler broadening and eliminates the dephasing by Wl ■ 
Thus, at the second window, the dressed atom-field dark 
state is stable with respect to the Doppler effect, and 
this stability enables observing the second window even 
at higher temperature. 


b. D s > D p 

For D s D p , the term VL v pn(t) — D*pi 4 (t) of 
Eq. (IA11I) is neglected because its effect is very small. 
Therefore, for this case, the coherence P 14 has negligi¬ 
ble effect on the atomic population of state |1). Conse¬ 
quently, the population of |1) is not affected by Doppler 
broadening at <5 P = 5 C . 

Under the additional constraint that <5 ps = 0, only ,034 
affects the population by reducing , 033 , and, as P 33 for the 
Doppler-free atomic-field system is almost zero, reduc¬ 
ing P 34 thus has no effect on the population of state |3). 
Hence, the population of each level (1A11I) at zero tem¬ 
perature will be the same as for the population at any 
higher temperature. At steady state, the atoms are all 
trapped in the state |'0d) for i5 pc = 0 and to |^>p) for 

S ps = 0. 


Appendix C: Absorption maxima and minima 

In this section, we show the calculations leading to 
closed-form expressions for the heights, or maxima, of 
the two absorption windows and also the nadirs, or min¬ 
ima, of these absorption windows. These expressions are 
derived first for the stationary atom and then generalized 
to the Doppler-broadened case. 
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1. Stationary atom 

Identifying the maximum height, Zi maXl and minimum 
height, /i m i rll of the z th window is subtle because the two 
Lorentzian transparency windows are cut asymmetrically 
into the overall Lorentzian absorption peak correspond¬ 
ing to zero-coupling field. First, we consider the i = 1 
case. 

The maximum is calculated by setting fl c = 0 = Q s 
and evaluating Eq. CD at S p = 6 C = 0: 


= p- 


^P 1 1 ~ P44 
74 


(Cl) 


The minimum h m ; ni is determined by setting Q s = 0 but 
^ 0 and evaluating Eq. ([Till at <5 p = <5 C = 0. We 
obtain 


l _ V (pu ~ P 44 ) 72 T2~^° pv 

^mini “ ' U 


7472 + | fie 


(C2) 


with zero absorption attained for 72 = 0. If 72 7 ^ 0 but 
condition (flSl) holds, 


V (Pn ~ P44) 72 

|^c| 2 


(C3) 


and minimum absorption is reached. 

The maximum and minimum of the first transparency 
window are used to calculate the half-maximum 


ki := 


Zlmax, + h m ^ ^ ~ (27472 + |» c | 2 ) 


274 ^7472 + |n c | 2 ) 


Applying condition (fT5l) yields 

P ( P 11 ~ P44) 


= 


274 


(C4) 


(C5) 


For the i = 2 case, we have a Lorentzian transparency 
window cut into the absorption curve corresponding to 
the two conditions fl s = 0 and fl c ^ 0 holding. For our 
case of DDEIT, we set fl c = 2 6 S , which establishes the 
second transparency window centered at <5 p = i5 s , which is 
the point that the maximum peak height Zi maX2 occurs. 
Therefore, Zi maX2 can be determined by calculating the 
fl s = 0 absorption curve value at S p = S s : 

(Pn — P44) (74 + 

Zlmax 2 = V ---fA (C6) 

7 4 + 4 {Wc -<J s) 

Linder the approximation that 

|f\| 2 72 


74 > 


261 2 ’ 


we obtain 


= p- 


P 11 — P44 


(C7) 


(C8) 


Thus, Zi maXl ss Zi maX2 . 

In order to calculate the minimum of the second trans¬ 
parency window, we set fl s ^ 0 and fl c ^ 0 and evaluate 
Imyp 1 ) from Eq. dill) at <5 P = i5 s to obtain the minimum 


Zimiri2 — P 


(/Oil — P44 ) 73 . {P44 - P33 ) |Us 


7473 + |^s| 2 r 43 (7473 + |L2 S | 2 ) 


(C9) 

The first term in the right-hand side of Eq. dC9l) repre¬ 
sents the absorption minimum, whereas the second term 
represents the maximum gain (negative absorption). 

If we wish to reduce absorption, decay from level |3) 
must be minimized, i.e., 73 — > 0. For the case 73 ^ 0, 
condition d!6D must be satisfied to minimize the absorp¬ 
tion. As r 43 ~ 74 , Eq. dC9l) is simplified to 


= P 


(pil — P44) 73 . P44 — P33 


10*1 74 j 

As P 44 = 0 is assumed, gain exists only when 
|P44 - P33I |pil - P441 73 


74 


|Hs 


or, equivalently, if 


In our case 


|p44 - P331 7374 

Iph P44I |fl s | 2 

\P44 - P33 


\pil - P44| 


1. 


(CIO) 


(Cll) 


(C12) 


(C13) 


Therefore, we also require condition Cl in order to 
achieve gain. The half-maximum is then 


x 2 =p 


(P11 ~ P44) ^27473 + |fi s | ^ 
74 (7473 + |0| 2 ) 


(C14) 


+ 


(P44 — P33) IO 


r 4 3 (7473 + IOI 2 ) 
For condition (ED and 74 73 , 

{Pn ~ P33) 


= P- 


74 


(CIS) 


For pn ss P 33 , x? is located at zero absorption. 


2. Doppler-broadened susceptibility 

In the case of Doppler-broadened susceptibility, the ab¬ 
sorption profile near the center (corresponding to zero 
velocity) is quite flat. As the two windows occur near 


74 
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the center; the maxima for both windows are the same. 
The maximum value is calculated for fi c = a = 0 and 
for (i p = S c : 


^Dmaxi 2 — 


7 ?a/ 7 T (pil ~ £44) 
74 + W L 


(C16) 


The minimum value of the first window is calculated for 
a = 0 and (5 p = ii c : 




- ^44)72 


72 “►O 


|fic 


■ 72 (74 + Wl) 


>0, 


(C17) 


which requires the condition |fi c | O 72(74 + !%) to 
hold in order to reach minimum absorption. 

For 72 yf 0, however, for Wl 74 this condition can 
be reduced to |a| 2 72 Wl- If the intensity of the 

driving field eliminates the inhomogeneous broadening 
due to Doppler broadening, it certainly eliminates the 
homogeneous broadening as well. The half-maximum of 
the first window is then equal to 


XTl = 


(pu - £44) 272 W L + |f2 c 


72 %l + 10:1 (74 + 1 %) 


(C18) 


For i = 2 , the minimum is calculated for fL = 0 and 


<5 P = 8 S with the result 


^Dmin2 — Tjy/ll 


(,£11 ~ £ 44 ) 73 
73 (74 + 1 %) + |0 
(£44 ~ £33) |Q|" 

73 Wl + (Wl - 74) IO 


2 (£44 — £33) |a| 2 

W L (27374+ IO I 2 ) 


(C19) 


The first term in the right-hand side of Eq. (1C 191) repre¬ 
sents the absorption minimum. This term tends to 0 if 
73 ~+ 0 . 

For the case of nonzero dephasing or relaxation decay 
from state |3), condition |fl s | 73 (74 + Wl) is required 
to minimize the absorption. For Wl 74, this condition 
can be reduced to m- The last two terms of the right- 
hand side of Eq. (IC19I) represent the maximum of the 
gain. After solving some algebraic expressions, Eq. (IC19I) 
becomes 


^]Jmin2 —£ 


(£11 — £44) 73 
73 (74 + Wl) + |fi s 


(C20) 


(£44 — £33) |a | 2 

Wl( 2 7 374 + |D s | 2 ) 

,, 2 73 W l (74 - W L ) + |D S | 2 ( 2 74 - W L ) 

73 Wl - |D S | 2 (74 - W L ) 

Now, we want to examine whether condition dill) , for 
73 0, is required to observe gain of the Doppler broad¬ 

ening susceptibility. If not, then the existence of gain 


suppresses absorption, and the second transparency win¬ 
dow is observed even if condition (1411) fails. We evaluate 
Eq. (IC20I) for the condition Wl » 74, in order to simplify 
the calculation, and evaluate for condition (1IB1) . which is 
necessary to minimize the absorption as shown earlier: 


I^Dmiri2 — 


(£11 — £ 44 ) 73 


73 W L + |fi, 
(£44 - £33) 


(C21) 


W L ^ 

In order for gain to exist, 

I£44 £33 1 I * , 73 Wl 


1 + 


73 W L 


73 W L + | a 


W L 


1 + 


7sW L + ia 


> 


|£n “ £441 73 


73 W L 


iaf 

(C22) 


which can be simplified by rearranging terms and substi¬ 
tuting the quantity 


to yield 


I £44 - £33! 

l£n “ £441 


iai 2 + w L 

73 Wl 


1 


» 0. 


(C23) 


(C24) 


Condition (1C24I) is always valid even if condition (l4lT) is 
not satisfied. Note that the derivation of inequality (IC24I) 
is based on the validity of condition (TTCJl) for homogeneous 
broadening. Therefore, condition (1161) is required for the 
gain to exist in our system, whereas condition dUD is not. 

The half-maximum for the second EIT window, with¬ 
out making any approximation, is 


>C2 =- 


(£11 - £44) (273 (74 + W L ) + |a| 2 ) 
(73 (74 + W L ) + |a| 2 ) (74 + W L ) 

(£44 — £33) iai 2 


W L (2 7 374 + |a| ) 


. 27 3 Wl(74 - Wl) + iai 2 (274 - Wl) 


-73W L 2 + iai (74 - Wl) 

Applying condition (fTTH) simplifies this expression to 

nrf ; (£n - £44) (273WL + iai 2 ) 


(C25) 


>r 2 =- 


( 73 WL + iai 2 ) (74 + Wl) 
^ 2 73 W l (74 - W L ) + |a| 2 (2 7 4 - W L ) ' 


-73 Wl + |a| (74 - W L ) 


. £44 ~ £33 
W L 


(C26) 
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which is the same as Eq. (l46l) for pn = p 33 « 0.5 and 

P 44 = 0 . 

For the condition that W\ J 74, Eq. (IC26I) reduces to 


^2 = 


vV^ ( 2 73 W L + |ft s | 2N \ ^ 

m(^7iaFj (wi - ro) ' (C7> 


Thus, the half-maximum of the second EIT window for 
high Doppler broadening depends on the population dif¬ 
ference between states |1) and |3). For equal population, 
the half maximum is always located at zero where ab¬ 
sorption vanishes. 
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